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Abstract 

We construct local solutions to 11-dimensional supergravity (or M-theory), 
which are invariant under the superalgebra D(2, 1; d; 0) © D (2, 1; d; 0) for all values 
of the parameter d. The BPS constraints are reduced to a single linear PDE on a 
complex function G. The PDE was solved in [1] modulo application of boundary 
and regularity conditions. The physical fields of the solutions are determined by d, 
a harmonic function h, and the complex function G. h and G are both functions 
on a 2-dimensional compact Riemannian manifold. The harmonic function h is 
freely chosen. We obtain the expressions for the metric and the field strength in 
terms of G, h, and d and show that these are indeed valid solutions of the Einstein, 
Maxwell, and Bianchi equations. Finally we give a construction of one parameter 
deformations of AdS? x S 4 and AdS^ x S 7 as a function of d. 
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1 Introduction 

We construct the exact local half-BPS flux solutions of 11-dimensional supergravity which 
are invariant under the superalgebra D(2, 1; d; 0) © D(2, 1; c'; 0). A noteworthy feature 
of these solutions is that the D(2, 1; d; 0) superalgebra is the unique simple superalgebra 
with a continuous parameter, d . The existence of this parameter makes the family of 
solutions particularly rich and opens the door to finding interpolating families of solutions 
as a function of d . 

We reduce the BPS constraints to a single linear PDE for a complex function G, 
2d w G — (G + G)d w \nh, where h is a freely chosen harmonic function. The differential 
equation is independent of d. As a result, given a harmonic function h, one can solve for 
G along the lines of [1], where the solutions for three special values of d are discussed. 
However, boundary and regularity conditions must still be satisfied and they do depend on 
d . The physical fields of our solutions, i.e. the metric and the field strength components 
are determined by the choice of d, the harmonic function h, and the complex function 
G. h and G are both functions on a 2-dimensional Riemannian manifold, d is a freely 
chosen real parameter which determines the fermionic generators but does not affect the 
bosonic generators of D(2, 1; d; 0) [5]. We check that the Einstein, Maxwell, and Bianchi 
equations are satisfied for every such choice of c', h and G. It is convenient to define 
a function H related to G by a fractional linear conformal transformation involving d . 
Although the PDE is linear for the function G, the physical fields are simpler in terms 
of H (see sec. 5) and perhaps even more importantly, the range of H as a function of 
w, w is independent of d whereas this is not the case for G, meaning that regularity and 
boundary conditions are easier to apply on H (see sec. 5.2). 

Some of the original interest in these solutions was motivated by seeking solutions cor- 
responding to intersecting M2 and M5 branes. A lot of progress was made in finding these 



1 



solutions in [1], where a general ansatz for the supergravity fields and supersymmetry pa- 
rameters was proposed. In particular, the bosonic subalgebra, 5*0(2, 2) x SO (4) x S0(4), 
which is independent of d, is naturally realized on an AdS^ x S* 3 x S* 3 fibration over 
a two-dimensional base space S. The reduction of the BPS equations to 2-dimensions 
was carried out for general d. However, explicit solutions to the reduced BPS equations 
were only found for the special values of d G {1,-1/2,2} and it was not obvious that 
non-trivial solutions existed for general values of d . 

For the values d G { — 1/2, —2}, the superalgebra D(2, 1; d; 0) © D(2, 1; d; 0) is simply 
OSp(4*\2) © OSp(4*\2) which is a subalgebra of OSp(8*\4). This implies that the cor- 
responding solutions admit solutions which asymptote to AdS-j x S 4 , including AdS 7 x 
S 4 itself [5]. The geometry AdS 7 x S* 4 is the near horizon geometry of M5 branes, 
which are conjectured to have a dual description in terms of a six-dimensional con- 
formal field theory (CFTg) [13]. Although the CFT@ is not yet completely known, 
there has been work on understanding the theory in various limits (see for example: 
[14, 15, 16, 17, 18, 19, 20, 21, 22]). 

In [23] (see also [24, 25, 26]), it was argued that the CFT 6 admits self-dual string 
operators, which are higher-dimensional analogues of Wilson lines in gauge theory. In 
particular, there should exist self-dual string operators which preserve half of the super- 
symmetries, corresponding to the supergroup OSp(4*\2)Q)OSp(4*\2). These solitons arise 
from considering M2 branes ending on M5 branes, much in the same way one obtains 
Wilson lines by considering fundamental strings ending on D-branes. In [8], the dual 
supergravity solutions were found, which described general configurations of arbitrary 
numbers of M2 branes ending on M5 branes in the near horizon limit of the M5 branes. 
Using the extended solutions presented in this paper, one might hope to generalize the 
solutions of [8], which would imply the existence of additional operators in the CFT 6 
theory. In particular, the self-dual strings should come with an additional parameter 
corresponding to the choice of d. 

For the value d = 1, the superalgebra D(2, 1; c'; 0) © D(2, 1; c'; 0) is OSp(4\2, R) © 
OSp(4\2,H) which is a subalgebra of OSp(8\4,TL). This implies that the corresponding 
family of solutions admits solutions which asymptote to AdS^ x S* 7 , including AdS± x S 7 
itself. The geometry AdS 4 x S 7 is the near horizon geometry of M2 branes. Progress 
in the dual CFT description of M2 branes was initially made in [45, 46, 47, 48] and is 
commonly referred to as BLG theory. In [28], a generalization was proposed in terms of a 
supersymmetric Chern-Simons theory, which allows for an arbitrary number of M2 branes 
and is commonly referred to as ABJM theory. The reduced superalgebra OSp(4\2,H) © 
OSp(4\2,H), corresponds to deformations of BLG or ABJM theory by the insertion of 
1 + 1-dimensional interfaces or defects, which preserve half of the supersymmetries. 

In [7], a specific solution in the d — 1 class was constructed which is dual to a 
deformation of ABJM by a dimension two operator. This solution is reminiscent of 
the supersymmetric Janus solutions of IIB supergravity [29, 30, 31], whose dual CFT 
description is discussed in [32, 33, 34]. Progress in constructing the aforementioned 
dimension two dual operator on the CFT side was made in [35] for the BLG theory, but 
an explicit construction of the interface operator in ABJM theory has not yet been given. 

In [36], it was shown that there are no other asymptotic AdS± x S 7 solutions, besides 
the Janus solution, with the assumption that the two-dimensional base space, S, has disc 
topology. The absence of such solutions could be interpreted to mean that ABJM theory 
does not allow for more general defect /interface theories other than Janus, however, there 
has been work in studying such deformations on the CFT side. M2 branes ending on 
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M5 branes from the M2 brane point of view was first studied in [37]. Recent progress in 
studying such defects and interfaces has been made in [38, 39, 40, 41, 42, 43, 44]. The 
existence of such objects in the CFT should imply the existence of dual gravitational 
solutions. There are two ways to get around the null results of [36]. First one may look 
for solutions where E has other topologies, such as an annulus topology. Second, such 
interface and defect deformations might require one to consider values of d away from 
d = 1. 

One of the principal motivations for finding solutions for general values of d is the 
possibility to find new families of supersymmetric solutions which can be viewed as defor- 
mations of the maximally symmetric solutions AdS 4 x S 7 and AdS 7 x S* 4 . The existence 
of such families can imply the existence of corresponding operators in the dual CFTs and 
give predictions for how the operators can deform the CFT. For example, the interpolation 
from d = 1, corresponding to AdS 4 x S 7 , to d G {0, oo} requires the decompactification 
of one of the 5' 3 s at the endpoints of the interpolation. This implies that the Kaluza-Klein 
scale vanishes at this value of d. In the dual ABJM theory this would correspond to the 
closing of a "dimension gap" , in the sense that one would expect towers of operators with 
continuous dimensions. We present two explicit examples of such families of interpolating 
solutions at the end of this paper, one family contains AdS 4 x S 7 and the other family 
contains AdS 7 x S* 4 . 

Perhaps the most interesting possibility is to look for a family of solutions which 
interpolates from AdS 4 x S 7 to AdS 7 x S* 4 . These would be dual to a family of CFTs 
which would interpolate as a function of d from the 3-dimensional ABJM theory to 
the 6-dimensional CFT@. We note that any such interpolation must pass through a 
decompactification limit. To support this idea, we note that in [49, 50, 51], there has 
been some progress on the CFT side in extracting some of the known CFT 6 data from 
massive ABJM theory [53, 52]. 

The organization of the paper is as follows. In section 2, we give the ansatz and review 
the results of [1] which we use. In section 3, we reduce the BPS equations, starting from 
expressions obtained in [1] to a single linear PDE. The expressions for the metric and 
the field strength in terms of h and H are obtained in section 4. We also show that the 
solutions solve the Bianchi identities, as well as the Maxwell and Einstein equations of 11- 
dimensional supergravity. In section 5, we give a summary of the solution, discuss some 
of the general features and give examples of one parameter deformations of AdS 7 x S* 4 
and AdS 4 x S 7 . 

2 Setup 

In this section we first give the D(2, 1; d; 0) © D(2, 1; d; 0) invariant ansatz for the 11- 
dimensional supergravity fields and then discuss the reduction of the corresponding 11- 
dimensional BPS equations to 2-dimensions. This section is a review of results from [1] 
and for a fuller understanding, we refer the reader to sections 2-4 of [1]. The reader 
familiar with [1] can skip this section. 

We will look for bosonic solutions to 11-dimensional supergravity which is defined by 
the following action 

s - 4 / dU *M R ~ rs F ^ pMNPQ ) -jkI CAFAF (21) 
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As usual, Kn is the 11-dimensional gravitational coupling, which is simply the 11-dimensional 
Newton's constant times a factor of 16n, g is the determinant of the metric and R is its 
Ricci scalar, F is the 4-form field-strength and C is its 3-form gauge potential, F = dC. 
The equations of motion are given in section 4.3. In the above action, the fermionic field, 
the gravitino, has already been set to zero. This is always a consistent choice. However, 
in order for the solution to preserve supersymmetry, the variation of the gravitino with 
respect to some supersymmetry parameter e must vanish. We will call such an e an 
11-dimensional Killing spinor. It is a 32 component Majorana spinor in 11-dimensions 
and the vanishing of the variation of the gravitino, called the BPS constraint, is given by 
the equation 

v M e + ^2 ( - r A f(r • f) + 3(r • F)r M y = o (2.2) 

where Vm is the covariant derivative with respect to the Levi-Civita connection for the 
metric. Tm are the usual gamma matrices satisfying the 11-dimensional Clifford algebra 
{T a , r 6 } = 2r] ab I 32 , where r] ab is the Lorentz metric in 11-dimensions and / 32 is the 32- 
dimensional identity matrix. T ■ F means the contraction of a rank 4 anti-symmetric 
tensor of T matrices with the field strength. In general, there is a supersymmetry for 
every linearly independent e which satisfies (2.2). 

We will look for bosonic solutions invariant under the 5*0(2,2) x SO (A) x SO {A) 
isometry, as this is the maximal bosonic subalgebra of D(2, 1; c'; 0) © D(2, 1; c'; 0) for any 
value of d . This leads to the following ansatz. The metric has the form 



ds 2 = fl ds 2 AdS3 + fl ds 2 s s + fl ds 2 sl + dsl (2.3) 



where ds 2 AdSs is the unit radius metric on 1+2 dimensional anti-de-Sitter space while ds 2 s3 
and ds 2 s3 are unit radius metrics on the 3-sphere. The metric component ds| is a metric 
on a 2-dimensional Riemann surface S. 

It is convenient to introduce a frame e A = dx M e M A , with A = 0, 1, • • • , 9, t|, where 
\ denotes the 10th spatial direction. We also introduce the notation e* 1 , e* 2 , and e* 3 for 
the frames on the unit radius spaces AdS 3 , Sf , and S| respectively. The 11-dimensional 
frames may then be written in terms of e %1 , e 12 , and e* 3 as follows 

e n =he 11 h = 0,1,2 

e t2 =f 2 e h 2 2 = 3,4,5 

e l "=he li ^3 = 6,7,8 
e a a = 9,1] (2.4) 

In terms of the frame components, the ansatz for the field strength is 

T? — n P 012a -L- n P 345a -L- n ^ 678a (9 ^ 

with the corresponding gauge potential given by 

C = M 012 + b 2 e 345 + b 3 e 678 (2.6) 

We use the shorthand notation e l i- k = e l A e- 7 A ... A e k . 

We introduce complex coordinates w, w on E, so that the metric takes the form 
<is| = 4p 2 \dw\ 2 . In terms of real coordinates x and y defined by w = x + iy, the metric is 
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given by ds| = Ap 2 (dx 2 + dy 2 ). The frames e a may then be expressed in terms of x and 
y as e 9 = 2pcfe and = 2pdy. We also introduce a complex frame e z on S defined by 

= (e 9 + ie*)/ 2 = p dw e~ z = (e 9 -ie*)/2 = p dw (2.7) 

in terms of which the metric becomes gLs| = 4e z e z . The field strength components can 
be written with frame or coordinate indices as 

11. . , . 

9iz = 9i9 - = ~9iw = 7T\9ix - Wiy) (2.8) 
p ip 

In general requiring the existence of non-vanishing e which satisfy the BPS equation 
(2.2) places restrictions on the supergravity fields (see [3, 4] for a comprehensive analy- 
sis). In the case we consider here, the supergroup D{2, 1; d; 0) © D(2, 1; c'; 0) has sixteen 
independent supersymmetries and so we will require the BPS equation to have sixteen 
linearly independent solutions. In this case, the BPS equation will completely constrain 
the metric factors /i, f 2 , fz, p (up to a sign), and the field strength components g± a , g 2a , 
93a 1 - 

To implement the D{2, 1; c'; 0) © -D(2, 1; c'; 0) symmetry, we first decompose the su- 
persymmetry parameter e as a tensor product of Killing spinors on the unit radius space 
AdSz, Sf and S 2 as 



= E 



x" M ®(r, M (2-9) 

■^m,V2,V3 are 8-component Killing spinors on the symmetric part of the space, i.e. they 
are tensor products of the 2-component Killing spinors on AdS$, S 2 , and S%. There are 
eight x T ' 1 > ? ? 2 ''? 3 , labeled by ^ = ±, i] 2 = ±, and ^ 3 = ±. The eight corresponding Cm,m,m 
are 4-component spinors. From the Majorana condition on e, the spinors can be shown 
to obey the following reality conditions 

( X m r = (h © o 2 © a 2 )x m3 {U, m , m T = (<x 3 © ^)C m , V2 , m (2.10) 

Due to symmetries of the reduced BPS equations, the eight Cm,m,m can & U De related to 
each other by the following operations. 

So : C -»■ i(I ® °" 3 ) C which has the effect of ^ ->■ (2.11) 

5,- : C -> s^cr 3 ' © C which has the effect of ^ ->■ -(-) <5 '' J ^ and g ia ->■ —(—) Si - j gia 

where j = 1,2,3, with si = «, and s 2 = s 3 = 1. So, S'i, S 2 , and S3 generate the 
full symmetry group {/, So, Si, SoSi, S 2 , SoS 2 , S3, So S3} which relates all eight Cm,m,m- 
Moreover, the four components of C+++ (and any other Cm,m,m) are reduced by the reality 
condition (2.10) to two independent complex components, a and j3 



,+++ 



/ a \ 

a 

V J 



(2.12) 



1 We will be able to completely determine all the physical fields i.e. metric factors and fields strength 
components in terms of functions h and H (see (4.9)) using the BPS constraints only. However, it will 
turn out that the Bianchi identities and the Einstein equations will add the constraint that Im(H) < 0. 
It is possible that this constraint could be found from the BPS equations only, by plugging the solutions 
back into these equations. We did not attempt this. 
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The BPS equation in the directions of the symmetric spaces reduces to six algebraic 
conditions on a and /3. Three of these are solved to yield expressions for the metric 
factors in terms of a and (3 

A = -(H 2 + |/3| 2 ) / 2 = --(M 2 -|/3| 2 ) h = -{aP-aP) (2.13) 

C\ C 2 C 3 

where ci, c 2 , and c 3 are real integration constants. These integration constants are directly 
related to the d parameter of the D(2, 1; d; 0) algebra by the relation d = c 2 /c 3 [5]. Two 
more of the algebraic conditions lead to equations relating the fluxes g lz , g 2z and g 3z 

g lz (a 2 - (3 2 ) = 2(ci - c 3 ) + 2ig 3z a/3 

g 2z (a 2 + (5 2 ) = 2(c 2 - c 3 ) + 2* fe a/3 (2.14) 

The last algebraic condition yields the constraint 

ci + c 2 + c 3 = (2.15) 

As a consequence of this constraint, there is only one independent parameter out of the 
three q. The constraint removes one degree of freedom and an overall rescaling of all 
three parameters can be absorbed into the definition of a and /3. We will often use the 
parameter c — (ci — c 2 )/c 3 as the remaining parameter. This parameter is, of course, 
uniquely related to the d in D(2, 1; c'; 0) by the formula c = —2d — 1. 

The BPS equation in the directions along £ is reduced to the following differential 
equations on a and (3. 

i 1 1 i 

D z a = +-u z a + Y^9uf3 - —g2 Z fi - J^93 Z u 

% 1 1 % 

D z (3 = +-u z (3 + ^9i z a + Y2 92za + i2 93zf3 

D z a = --u z a + -g lz p - -g 2z p + -g 3z a 

i - 1 1 % 

D z f3 = - ^UzP + -^gizOL + -g 2z a - -g 3z (3 (2.16) 

where Cj z is the spin connection on E and we have introduced the notation 

D z = 9 - D,= ^ (2.17) 
P P 

as implied by (2.7). For vanishing torsion, the frame is related to the spin connection by 
e a + u a b A e b = 0, which can be used to compute Co z (the ^-component of w 8 9 ) in terms of 

P- 

.d w p A .dap 

u z = % — —aw oo z = —1 — —aw (2-lo) 
p 2 p 2 

Finally, it was shown in [1] that the system of equations (2.16) after using the algebraic 
constraint (2.14) admit a first integral for any value of c. We will put off giving this first 
integral until the next section where we will rederive it. In [1], the full solutions of the 
system of equations defined by (2.14) and (2.16) were found only for the special values 
of c = 0, 3, —3. The values c = 0, 3 correspond to solutions asymptotic to AdS-j x S* 4 , 
while the value c = —3 corresponds to solutions asymptotic to AdS^ x S 7 . However, we 
note that the reductions reviewed in this section and the form of the first integral are all 
valid for general c. In the remainder of this paper, we solve the system of equations for 
general values of c. 
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3 Solution for General Values of c 



In this section, we derive the solutions to the system of equations defined by (2.14) and 
(2.16) for general values of c. The generalization of the solutions of [1] to general c will 
require a modification of the methods of section 5-7 of [1]. 



3.1 Reduction to two equations and the first integral 

We begin by introducing a s and (3 S by rescaling the variables a and (3 in the following 

way 



1 //t 3 \ 8 1 f K 3 



a = — — a s a = — — a 
pz \ k J pi \ K 



P 4 \ K J P z 



(3.1) 



where k is an arbitrary holomorphic function, whose significance will be seen shortly 2 . 

Next we rewrite the differential equations (2.16) as logarithmic derivatives for the 
combinations oo a = a s a s 3 , oop = f3 s /3 s 3 , ou a = a s a 3 and Qp = f3 s /3 3 . Although, we could 
express a, (5 in terms of these new variables, we do not do so for now. These definitions 
are made by demanding that both d w p and g 3z drop out of the equations for u a , and up. 

d w \nu a = -p{gi z - g 2 z) V — 

4 \a s a s 

1 / \ f a s &s\ 

d w Inup = -p{gi z + g 2z ) I -j- + 5- 

I ~ (3 ' 2) 

d w \nu a = 2d w \np- d w \nn+ -(g lz - g 2z )p ( 3— + - — J +i-pg 3z 



a, 3 a,/ 3 



d w \nujp = 2d w \np- d w \nn+ ^(gn+g 2z )p (^j^ + \ ~ ^P93z 

We now proceed to exclude the g^s which will leave us with three equations. First, 
we rewrite the algebraic constraints (2.14) to have one equation without g 3z and another 
to be homogenous in the gi z (i.e. without the constant term) 

(9iz - 92 Z )a 2 ~ (9iz + 92 Z )P 2 = 2(ci - c 2 ) ^ 3 ^ 

g lz (a 2 - (3 2 )(c 2 - c 3 ) - g 2z [o? + (3 2 )(d - c 3 ) + 2ig 3z a(3(c 1 - c 2 ) = 

Next we obtain an inhomogenous differential equation without d w p and a homogenous 
differential equation. To do so, we use the first two equations of (3.2) to express the 
combinations (g± z + g 2z ) and (gi z — g 2z ) in terms of d w lnu a and d w Inup. To obtain the 
inhomogenous differential equation, we insert these expressions into the first equation of 
(3.3). To obtain the homogenous differential equation, we add the last two equations of 
(3.2) which eliminates g 3z and then eliminate the expressions for (gi z + g 2z ) and (gi z — g 2z ) 

2 This rescaling is inspired by the form of the first integral, known from [1]. However, we stress that 
in this equation k is not a first integral but simply an arbitrary holomorphic function. 
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(using the first two equations of (3.2) as before). The resulting two differential equations 
are 

3 

K ( p \ 2 - 

d w O a - d w Qp =-( Cl - c 2 ) —r (a s /3 s + a s /3 s ) 

1 VN/ _ _ (3_ 4 ) 

d w u a d w U(3 p 1 d w uj a 9/3 s a s + (3 s a s 1 d w u>p 9a s (3 s + a s (3 s 

+ =4d w ln— + - — — — =- + -— — — =- 

u a up \k\ 3 co a a s p s + a s p s 3 up a s p s + a s p s 

There then remain two more independent equations. One of them must give the 
expression for g 3z . g 3z is obtained using the difference of the u a and ujp equations, i.e. 
the last two equations of (3.2), but this equation is not displayed. We get the final 
equation by inserting this expression for g 3z into the second equation of (3.3) and once 
again insert the expressions for {g lz + g 2z ) and {g\ z — g2 Z ). The final result is a total 
derivative 

dj(a s a 3 s - M!) ~ C1 ^a s p- S {a s p- S - a A)) = (3.5) 

The above equation implies that the expression in parenthesis is a first integral. Using 
the fact d w R = 0, the above equation implies that 

<9„, jp((aa 3 - (3l3 3 ) - Cl ~ C2 a/3(a/3 - apf) } = (3.6) 

In other words, the expression in parenthesis above is an anti-holomorphic function. Now, 
since the function R used in the rescaling (3.1) was an arbitrary anti-holomorphic function, 
we can choose 

R = p{{aa 3 - P(3 3 ) - ( ^^a(3{ap- - aft)) (3.7) 

C3 

After this identification, we see that R is a first integral. This is also the expression for 
the first integral which was derived in [1]. Note that, as was observed there, k is a (1,0) 
form. Part of the convenience of the rescaling 3.1 is that while a and f3 are (—1/4, 1/4) 
forms, the quantities a s and (3 S , as well as their complex conjugates, are (0, 0) forms. 

Before we proceed, we would like to discuss how the case of c\ = c 2 (i.e. c = 0) fits 
into the story. The careful reader may have observed that in this case the inhomogeneous 
differential equation, (3.4), reduces to (3.5)! This is a consequence of the fact that the 
algebraic constraints (3.3) become degenerate. However, this is not actually a problem. 
Since this is simply an intermediate step in the calculation, we can proceed for general c 
and take the limit c — > once the full solution is found and recover the c = solution 
found in [1]. Indeed, it will turn out that all of the physical variables of the problem 
(metric factors, field strength components, spinor components) are expressed in terms 
of a harmonic function, h, and another function, G, which satisfies a simple partial 
differential equation. All of these expressions and the partial differential equation will be 
valid for any value of c. 

We note that it is also certainly possible to follow a derivation which is valid for any 
c at every step. A different inhomogeneous equation can be found by going back to an 
earlier form of the algebraic constraints (2.14) which are not degenerate. We do not 
pursue this course here because the inhomogenous equation obtained this way is more 
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complicated for general c and we already have a derivation for c = valid at every step 
in[l]. 

We now continue working with our equations. We can rewrite the constraint (3.7) in 
terms of the rescaled variables a s and (3 S , or alternatively in terms of uo a and uip 



1 = [a s a] - Ml) - 



ci - c 2 



c 3 



1 + 



gi - c 2 

c 3 



ci - c 2 

c 3 



W/3 1 + 



ujp 



(3.8) 



Since we have two constraints, (3.7) and its complex conjugate, we should be able to 
rewrite our system (3.4) in terms of only two variables. To this end, we introduce a new 
variable A 4 = w^/w a and its complex conjugate, in terms of which the second equation 
of (3.8) can be written as 



;i + cAA)-A 4 (l + c(AA)- 1 ) 



(3.9) 



where as mentioned before c = (ci — c 2 )/c3. This will allow us to eliminate ujp, u a , u a 
and up in terms of A and A. We can now easily express equations (3.4) in terms of u>p, 
u a ,u a , up using 



a s /3 s 



up 
up 



(3.10) 



The equations (3.4) in terms of A and u a and their conjugates are 



d w u a - d w (\ A uj a ) = |(ci - c 2 ) 
^-d w \n(\ 2 u a ) = d w \n ( ^~ 

4 \ AC 



2 i \- 

WoWa 



(AA)2 

1 , A 2 + 9A 2 
8^ ln ^A^P" 



1 



(3.11) 
9A 2 + A 2 



+ -d w ln(A Q a ) A2 + A2 



This system of equations can now be expressed in terms of only A, A, and p by plugging 
in (3.9) and its complex conjugate into (3.11). However, we do not write down these 
expressions explicitly because they are unwieldy and unilluminating. Instead we now 
explain in general terms, the strategy we use to perform an integration, executing the 
analytic but messy algebraic computations in Mathematica. 



3.2 Integrating out p 

The system of equations (3.11) can be put into the following general form 



C(\,\)d w \ + D(\,\)d w \ 



ci - c 2 p 
2 V « 



A(X, X)d w X + B(X, X)d w X = d w In 



K\ 



(3.12) 
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We take inspiration from section 5 of [1] where we see that p is multiplied by some 
multiplier M(A, A) so that it is possible to integrate out a combination p and M by 
equating the differential forms on the left hand sides of the resulting equations, and 
thereby, also equating the right hand sides of the resulting equations, which contain 
only p and M with no explicit A, A dependance. In other words, we multiply the first 
equation by a function M(A, A), while in the second line we add to both sides the quantity 
d w lnM(A, A) forming the following equivalent system 



M(Cd w X + Dd w X) = %r^n (ri) M 

2 \ K 



(A + d x \n(M))d w X + {B + d- x \n(M))d w X = d w In (m (j^j ' ) (3.13) 

We now demand that the left hand sides are equal and moreover that the coefficients of 
d w \ and d w X are equal. If this demand is satisfiable, we can equate the right hand sides 
as well and get the following easily integrable equation. 



M 



3 



Note that there is no guarantee that an M can be found which would equate the left 
hand sides of (3.13) so we will show that in this case, such an M indeed exists and then 
we will return to integrating (3.14). Equating the coefficients of d w X and d w X in (3.13) 
yields 

19M A = CM L d JL +B = DM (3.15) 



M dX M OX 

Fortunately, we do not even have to solve these differential equations to find M (assuming 
it exists) since it must satisfy both differential equations and is therefore overdetermined. 
We can show that if M does exist, it can be solved for algebraically. Although, then it 
must still be tested whether (3.15) are solved. To derive the algebraic expressions, we 
linearize (3.15) by dividing by M and rewriting the equation in terms of the new variable 
N = 1/M 

ON dN 

+ AN = C — + BN = D (3.16) 



dX dX 

We then differentiate the resulting equations such that both contain the second order 
mixed derivative of N, and use (3.16) to exclude the first order derivatives of N. The 
second order derivative cancels between the two equations and the following expression 
is obtained 



1_ AD-BC + ^g-^ 
M ~ 



d _A_al 9X (3-17) 
ax ax 



Applying this procedure to (3.11) with (3.9) inserted yields the following expression for 
M 

i _ c (a 2 -A 2 )(a 2 A 2 -i) f xx y 

(3.18) 



M 2 (AA)I V(A + cA 2 A-cA 3 -AA 4 )(A + cA 2 A-cA 3 -AA 4 ) 
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We note that M is imaginary. This ugly expression will soon be shown to be something 
quite nice, but first, we discuss what the equations have become. The equation (3.14) for 
p can now be easily solved by rewriting it as 



(3.19) 



Since k is holomorphic, we can integrate both sides. To do so explicitly, we introduce the 
imaginary harmonic function h by 




d w h = — k (3.20) 



Integrating, we then have 



\K\ 



3\ -1 



= (3.21) 



Finally we are left with one very ugly equation equation to solve. To obtain it, we take 
the second line of (3.12) and use (3.21) to eliminate p 

d w \nh = A 2 d w X + B 2 d w X (3.22) 
111 1 2cA(A 2 + A 2 ) 

^2 = -T T + T T + TT7T 77 ~ 777T 77 + 



A - A A + A A(AA-l) A(AA + 1) A(-A - cA 2 A + cA 3 + AA 4 ) 
11 A A 2(cA 2 - 3cA 2 - 4AA 3 ) 4A 

E>2 = t r + T r + TT 7 + TT 7 H ; 777 7^ 7T7 + 



A-A A + A AA-1 AA + 1 -A - cA 2 A + cA 3 + AA 4 A 2 + A 2 

Since k is an arbitrary holomorphic function, h is an arbitrary harmonic function that 
we choose. The choice of h determines k and the above equation is solved for a given h. 
This final equation may then be viewed as a differential equation for A and A. 



3.3 Linear equation 

Our next goal is to integrate (3.22). To start, we map it to a linear equation. The general 
idea, which was carried out for the special values of c = —3,0,3 in [1] is as follows. Life 
would be easy if the right hand side of (3.22) would be a total derivative, but this is not 
the case. The next best thing that can be tried is to multiply (3.22) by some multiplier, 
m, such that the right hand side with this new multiplier is a total derivative of some 
function, which we call G(X, A). To find such a multiplier one would have to solve the 
equation 

d(mA 2 ) _ d(mB 2 ) 

~ dX " dX {6 - 26) 

Not only is this quite difficult, but in and of itself, would not directly lead to the integra- 
tion of (3.22). Instead we look for a multiplier which linearizes the equation. This means 
that m should be a linear function of G and G. This lucky state of affairs was indeed the 
case in [1], where the expressions A 2 and B 2 were much simpler and m — (1/2) (G + G). 
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We try our luck and assume for general c, that m also has the form m — (1/2) (G + G) and 
we can obtain the same integrable equation which was obtained in [1] for the G function 3 

d w G= \{G + G)d w \nh (3.24) 

We describe our procedure to integrate the equation in general terms but the algebra 
is performed in Mathematica. Before we start, we note that, as can be seen below in 
(3.27), the numerical coefficient in front of G + G drops out of the integrability conditions 
and one can look for m in the more general form m = constant * (G + G). However, the 
procedure picks out a unique coefficient when the final expression for G(A, A) is plugged 
into (3.24) and compared with (3.22). For our case, one finds that the coefficient is indeed 
1/2. After multiplying (3.22) by (1/2) (G + G) and equating the result with (3.24), one 
obtains 

d w G =d x Gd w \ + d x Gd w \ = ^A 2 (G + G)d w \ + ^B 2 (G + G)d w \ (3.25) 

Equating the coefficients of d w X and d w \ we obtain the conditions 

d x G = ^(G + G)A 2 d- x G=±(G + G)B 2 (3.26) 

The system of equations, (3.26), is an overdetermined system. Similar to the case of 
looking for M, once we find G, using the integrability conditions for the system (3.26), 
we must check that 2d w G/(G + G) is actually equal to the right hand side of (3.22). To 
derive the integrability conditions, we take the mixed derivative of (3.26) and equate the 
second order derivatives. We also note that the first order derivatives are simply related. 

dxd- x G = ^(A 2 l -(G + G)) = §- x (B 2 \(G + G)) (3.27) 

dxG = ^d- x G (3.28) 

D 2 

Note also that both equations can be complex conjugated providing additional equations. 
Solving these four equations together, we obtain 

d x G = (G + G) B ff 2 J*£ ] = (G + G)G 2 

d x G = (G + G) A ff 2 -°f 2) = (G + G)D 2 (3.29) 

A 2 A 2 — D 2 ti 2 

where we have defined new variables C 2 and D 2 . Adding the above equations, we obtain 
a simple equation and solution for G + G. 

|-(G + G) = (G + G)(G 2 + J D 2 ) 

OA 

G + a = »(A«-A«)(1-A*A°) 

" 1 (A(l-A 4 ) + cA(A 2 -A 2 ))(A(l-A 4 ) + cA(A 2 -A 2 )) 1 ' ; 

3 Note that more generally one could try the ansatz m = aG + bG with a and b arbitrary complex 
constants. However, this is more complicated and not necessary for the problem considered here. 
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where /(A) is an arbitrary function which has appeared as a result of the integration. 
Fortunately, we are able to remove this degree of freedom because the expression for 
G + G has to be real. The result of the integration actually has some branch cuts and the 
correct choice using this function is just to have the % in the numerator as shown. We still 
have the freedom to multiply by an arbitrary real constant but that is just a symmetry 
of equation (3.24). To find G — G we must subtract the equations of (3.29) and integrate 
again. 

3 (G-G) = (D 2 -C 2 )(G + G) 



r r _ -*(A 2 -A 2 ) 2 (l + 2cAA + A 2 A 2 ) 

(A(l-A 4 ) + cA(A 2 -A 2 ))(A(l-A 4 ) + cA(A 2 -A 2 )) 1 ' ' 



Note that we had the freedom of an additive function of A in the second equation of (3.31) 
coming from the integration, but as before this freedom is reduced to that of adding a 
constant because we know that this expression has to be pure imaginary. This constant is 
fixed to be i/c and is already included in the above expression. We fixed this constant by 
demanding that the expression has a finite limit as c — > since we know the solution to 
be perfectly well behaved at this point from the results of [1]. We also note that equation 
(3.24) allows us to add an arbitrary imaginary constant to G. G is given by the following. 

(A>A)A(A-A)_ 
A + cA 2 A - cA 3 - AA 4 K J 

With this identification we can check that (3.32) inserted into (3.24) exactly reproduces 
(3.22). Since this equation is independent of c and identical to what was obtained in [1], 
the integration method used there (in section 8) to integrate (3.24) applies to our work 
as well. However, the range of G and the expressions of the physical fields in terms of G 
do depend on c as is discussed in section 5.2. 



3.4 Comparison with previous results and discussion 

We would now like to compare our results to those of [1] for the special values of c that 
are solved there. To this end, we will express G in terms of the original variables a, (3, p. 

ia/3(a/3 — a/3) 

G = — ^— — — =r = ia/3{afi - a/3)p/R = ia s /3 s (a s /3 s - a s /3 s ) (3.33) 

aa A — p/3 6 — cap(pa — a/3) 

where to obtain the second equality we notice that the denominator is R/p. We must also 
express the G obtained for c = and c = — 3 in terms of a, (3, p. In the case of c = we 
get a straight-forward agreement. For the case of c = —3, we find that G appearing in 
[1] has the following form as a function of a, (3. 

. da 3 — (3(3 3 + a/3(a/3 — a/3) .(aa 3 — (3(3 3 + a/3{a/3 — a/3)) . . 
c= - 3 ' [1] = l aa 3 -(3l3 3 + ^a(3{(3a-a(3) = % W~P ' 

This expression looks quite different from (3.33) but it is actually equivalent up to a real 
multiplicative constant and an imaginary additive constant. As we have already noted, 
adding an imaginary constant or multiplying by a real constant is a symmetry of G since 
solutions differing in this way are not distinguished by (3.24). (We also note here that 
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due to branch cuts we can only claim to be able to express G in terms of a, (3 up to 
a multiplicative 4th root of 1. This is also true for the G expression (3.33). However, 
once we obtain the metric factor formulas, we will see that the above choice of branch 
cut yields AdSj x S 4 and AdS± x S 7 .) To show the equivalence of (3.34) and (3.33) we 
rewrite (3.34) by adding and subtracting i. 

_ . (aa 3 - fi(3 3 + a/3(aP - a/3)) - (aa 3 - fi(3 3 + 3a(3(f3a - q/3)) . 
c= ~ 3 ' [1] ~ 1 aa 3 - (3(3* + 3a(3{(3a - a(3) + % 

= -Aia(3((3a - a(3)p/R + i = -AG + i (3.35) 

We would now like to emphasize that amazingly, not only is G of the form (3.33), which 
is rather simple even for general c, but also the dependence on c can be hidden inside of 
p or inside of the scaled variables a s , (3 S . It is also interesting to look at the form of M 
in terms of a and (3. 

!_ = c(aa + Pf3)(aa - 0)(af3 - aft) 



2^(aa 3 - (3(3* - ca(3(a(3 - a(3))(a& - (3(3* - ca(3(a(3 - a(3))^j 

2(|«|/p)i 



(3.37) 



where to get the last equality we have used the expressions for the metric factors (2.13). 
Again, it is remarkable how simple this expression is and again the dependence on c can 
be hidden inside of p. Now if we remember the relationship between h, M, and p, (3.21), 
and we insert (3.37) we uncover another constraint. 



2 1 (\k\ 

h 



ci - c 2 M \ p 

iccic 2 c 3 



ci - c 2 



/1/2/3 = icic 2 fif 2 f 3 (3.38) 



Note that, as was already discovered in [1], the product of the metric factors is a harmonic 
function. Also note that we can rediscover the definition of k and get an independent 
check of (3.20). If we differentiate h using the second equality of (3.38) and use only the 
original equations (2.14-2.16), we recover k. 

d w h = dvC^^-f&h) = -P(aa 3 - (3P 3 - caP(aP - a(3)) = -« (3.39) 
c i _ c 2 



4 Derivation of Physical Fields 

In this section, we work out the expressions for the metric factors and fluxes in terms of 
the quantities G and h. 



4.1 Metric factors 

Our first goal is to find the expressions for the metric factors in terms of G and h. It turns 
out to be technically difficult to directly express a, (3 in terms of G, G so once again, we 
take inspiration from [1] and we express a, (3 in terms of trigonometric functions of if, (p 
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which can be expressed in terms of G, G. We define a new variable k which is simply 
related to k 

k = p(aa 3 - P/3 3 ) (4.1) 
R = k — cpa(3(a/3 — aft) 

Furthermore, consider the algebraic constraint given in the first equation of (3.8). The 
second term a s (3 s (a s (3 s — a s (3 s ) is just iG. We then introduce F and rewrite (3.8) as 
1 = F — (cG/i) where 



F = 1 - icG = a s a 3 s - /3 S (3 3 S 

k = RF = R(l - icG) (4.2) 

Note that F and k are just 1 and k for the special case of c = 0. In complete analogy 
with [1] we define <p, <p from the first equation in (4.1) as follows. 

sinh%) = ^(3(3 3 cosh 2 (v?) = {ao? (4.3) 

The complex conjugate expressions are also implied as usual. The convenience of tp is 
that these expressions are easily invertible and can be solved to obtain the following 

Now we express G in terms of if starting from (3.33) and using the above equations 
to eliminate a and (3. The resulting expression for G is 



From the above and its conjugate we obtain 



-GG 
FF 

H _GF _ sinh(2<^) 
H ~ GF~ sinh(2v?) 



V^HB = \l = sinh(v, - <p) (4.6) 

(4.7) 



where we have defined H = G/F. For completeness, we also provide the following easily 
derivable but useful formulas 



cosh((^ — if) = a/1 — \H\ 2 
sinh(2(v? - (p)) = 2 V /\H\ 2 {\H\ 2 - 1) (4.8) 

Note that, since f — (p is pure imaginary then | cosh(</? — <p)\ < 1 which implies \H\ < 1. 
Since G and H are related by H = G/F with F given by (4.2), or equivalently 

h -T^g < 49 » 

the restriction on the range of H implies a restriction on the range of G. This is further 
discussed in section 5.2. 
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We will need the expression for sinh(2</?), which requires quite a bit of algebra using 
the expressions (4.6) and (4.7). 

sinh 2 (2y?) = 4 ^ 2 (l g J 4 -| g | 2 ) (4 . 10 ) 

Finally, in the derivation of the currents, we will need sinh(</? + <p>) which can be obtained 
from (4.6), (4.8) and (4.10). Once again, branch cut choices are a difficulty and the sign 
between the two roots which appears in the expression below is chosen to be negative 
since sinh(<p + (p) vanishes for imaginary p or equivalently imaginary H. 

sinh(v? + 0) = \J-\H\ 2 {1 + sinh 2 (2^)) - ^ (1 - \H\ 2 ) sinh 2 (2<£) 
\H\(H + H) 



^-{H - H) 2 -A\H\^ 



(4.11) 



We can now obtain p from the first line of (3.38), (3.37), and (4.4). Note that in 
(3.37), the factor of (|/t|/p) 3 / 2 in the denominator is exactly what is needed to write M 
entirely in terms of p and p using (4.4). Then p is obtained by plugging M(ip,<p>) into 
the first line of (3.38). 



r — I — 1 

v2c 3 /i | sinh(2<^) 1 2 
Using (4.2, 4.8, 4.10) to eliminate p in terms of H we obtain 

« kl 6 |F| 6 



fc|s sinh(2(v? - y?)) 



P = 



4h* 



(\H\ 2 - l)((H - H) 2 + A\H\^ (4.13) 



ip = p — p e = 



Note that the first parenthesis is always negative since as we already argued \H\ < I. 
The second parenthesis is also always negative. To see the latter note that we can follow 
the logic of [1] once again and say that 

(H - H) 2 + 4\H\ 4 = -4sin 2 (p) sin 2 9 + 4sin 4 (/i) 

'sinh(2£>)y 
^sinh(2</?) J 

From the above definitions of 9 and p in terms of p and p it can be shown that 
tan#tanh(<p + <p>) — tan p. (This formula is also in [1] since the relation between 9, 
p, and p is unchanged. ) Using also that | tanh(</? + (p)\ < 1, the former implies that 
I tan// 1 < I tan 6*| which is equivalent to |sin/t| < |sin#|, implying that (H — H) 2 +A\H\ A < 
0. 

Now we find f\ by plugging (4.4) into (2.13) 

1 / I Fk\ \ 5 

h = -(\a\ 2 + \(3\ 2 )= V (Icosh^l + |sinh^|) (4.14) 



ci" ■ \ c\p 



By using the following trig identities 



cosh(/?| + I sinh</?| = ^| sinh(2</?)| ± y | sinh(2(/?)| 2 + cosh (p — p) j (4-15) 
cosh(/?| — I sinh(/?| = ^— | sinh(2</?)| ± \J \ sinh(2(/?) \ 2 + cosh 2 ((/? — p)^j (4.16) 
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and plugging in (4.8, 4.10) we obtain 



| cosM + I sinh v?| = ( {H ^2 + 1] m4 ) 4 (2|# I' ± KH ~ H)) h (4-17) 

| cosh <p\ - | sinh <p\ = ( (^jf^^l ) 4 (-2|#| 2 ± " ( 4 - 18 ) 
Finally plugging (4.13, 4.17) into (4.14) we get 

fi c\ 9/l „ . (2\H\ 2 ±i(H - H)) 3 

The denominator is always positive. The first parenthesis in the numerator is always neg- 
ative as has been already discussed and h 2 is also always negative since h is an imaginary 
function. The sign in the second parenthesis depends on whether the imaginary part of 
H is positive or negative and must be chosen such that the paranthesis overall is positive, 
in other words, such that the two terms in the parenthesis add. The reason for this is 
that this choice of sign exactly picks out | cosh tp\ + | sinh tp\ over | cosh tp\ — | sinh (p\ which 
picks out fi as opposed to / 2 . Also note that although there are no roots appearing in the 
expression above, square roots appear during the derivation so there is an overall choice 
of sign in the above expression. The positive sign is chosen to give a positive expression. 
Similarly, we also have 

h = — (M 2 -|/3| 2 ) = - V (|cosh^| - |sinh^|) (4.20) 

C 2 V C 2P J 

Plugging in (4.13, 4.17) we get 



c 



2 



(2! 


H 








-H) 2 + A\H, 


4\2 



H = ~4h\\H\ 2 - 1) ' ' '\ ' , (4.21) 



This expression is very similar to above, but here the sign choice in the second paren- 
thesis of the numerator is such that the two terms subtract, which is the sign choice 
corresponding to | cosh — |sinh</?|. As we have argued above |Im(if)| > \H\ 2 so this 
parenthesis is negative and the overall expression is positive. 
Finally we come to the last metric factor 

h = —(a0-M = -(—)* (4.22) 



Cs c 3 V P / sinh(2v9)|5 

Plugging in (4.6, 4.10, 4.13) we get 



h 2 {{H - H) 2 + A\H \ 4 ) 

4 (l#l 2 -i) 2 



Note that the sign choices inside the numerators of f\ and / 2 are correlated not only 
within the numerator but also between f\ and / 2 themselves. The sign is chosen such 
that ±i(H — H) = \H — H\. This implies that we can get rid of the sign choice in the 
above expressions which are in terms of H and write them with a unique sign using the 
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absolute value function. In the next section, we shall see that requiring the fluxes to 
satisfy the Bianchi identity selects the branch with i(H — H) = \H — H\, which means 
that lm(H) < 0. 

Taking c = and G = H in the final expressions for the metric factors (4.13), (4.19), 
(4.21) and (4.23), we recover the expressions for the metric factors obtained in section 5 
of [1]. Furthermore, if we take c = — 3 and G — > — (1/4)G + (i/4), we recover the metric 
factor expressions in section 7 of [1]. As discussed in detail in section 3.4, one needs to 
perform a linear transformation on G in order to take into account the difference in G 
functions in [1] and our paper. 



4.2 Field strength components 

Our next goal is to compute the field strength components, gi z . We will first compute 
gi z and g 2z . To do so, we use the first two equations of (3.2) to express g\ z and g 2z in 
terms of Cu a , Cup, a s , (3 S and their complex conjugates. Next we use the definitions of 
u a and up given above (3.2) and then (4.3) to express u a and Cup in terms of F and tp. 
Finally, we eliminate the remaining a s and (3 S , by first rewriting them in terms of a and 
(3 using (3.1), then the relations between a, (5 and tp given in (4.4) and finally we use 
(4.2) to write k in terms of F, with the overall k dependence now dropping out. Putting 
everything together, we obtain 

= 2 r vy~~~rr (d w log(Fcosh 2 tp)\ cosh tp\ + d w log(F sinh 2 tp) | sinh <p\) 

smh^ + ^ i (424) 

2 ^ smh ( 2 09 ) ^ 

P92z = . w — — ^(-<9 w log(Fcosh 2 (^)| cosh^l + d w log(F sinh 2 (p)\ sinh p\) 
smh((/9 + ip) 

The various quantities appearing in the above formula are given as follows: | cosh<^| and 
| sinh can be computed using (4.17), sinh(</? + tp) using (4.11) and |sinh(2</?)| using 
(4.10). To get cosh 2 tp and sinh 2 tp, we use the double angle formulas to express them in 
terms of cosh(2<^), which in turn can be expressed as a root of l + sinh 2 (2<^) with sinh(2<^) 
given by (4.10). 

We continue the calculation by switching to computing d w bi = —pffg iz /2, which are 
the quantities appearing directly in the Bianchi identities and equations of motion. 4 We 
make use of equations (4.19) and (4.21) for the metric factors. It is also necessary to 
remember that F can be expressed in terms of H by F = (1 + icH)^ 1 . In deriving the 
final formulas, the following algebraic identity is frequently used 

(2HH ± i(H - Sfj (2HH T i(H - Sfj = AH 2 H 2 + (H - H) 2 (4.25) 

It is important to note that in the expressions for cosh 2 tp and sinh 2 tp, a branch cut choice 
must be made and can be be determined by demanding that the Bianchi identities are 
satisfied. This is also the sign choice for which the expressions reduce to the ones in [1] 
for c = 0. Parameterizing, the branch cut choice for d w b\ by introducing v — ±1, we 



4 The factor of a half comes from our conventions for which g z = g 8 — igg and d w = (d x — id y )/2. 
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have 



d w h =ic 3 h(l - HH) ^cl(H + H) {{H - H) 2 + AH 2 H 2 ) 



Acd w H{2HH + iu(H - H)f 
-i + cH 



+ (hH{HH -1)(H-H + 2H 2 H)^j 1 (d w HH(H + H) + d w HH(H - 3H + AHH 2 )} 

- 2iV{H -H + 2H 2 H)(2HH + iu(H - H)) 
2(H -H + 2H 2 H) 2 (2HH + %v{H - H)) 2 ~ 



+ 



{H - H) 2 + AH 2 H 2 



(4.26) 



We have also introduced v = ±1, which is given by the sign of Im(— H) and parameterizes 
the sign choice in the metric factors. Demanding the Binachi identity to hold, one finds 
that v = —v. The same choice arises in the calculation of d w b 2 and the sign must be 
chosen in the same way. There is another sign ambiguity, since ff and /f are determined 
only up to a sign, we parameterize these sign choices by including overall factors v\ = ±1 
and i>2 = ±1. The final expressions for d w b\ and <9 W &2 are given by 



d w h =v x 2c z h c\{-i + cH)(H + H) [2HH - vi{H - H) 



Hd w H(-i + cH)(H -3H + AHH 2 ) 



- iud w H^ - iHH{H - 3H + AHH 2 ) + c((H - H) 2 - HH 3 + 3H 2 H 2 ^ 



+ Hd w H{-i + cH)(H + H)(1- ivH) 



(4.27) 



d w b 2 = - v 2 v2c z h \ c\(-i 



+ cH) (H + H) (2HH + ui(H-H)^j ^ 



Hd w H(-i + cH)(H -3H + AHH 2 ) 



+ iud w H^ - iHH{H - 3H + AHH 2 ) + c((H - H) 2 - HH 3 + 3H 2 H 2 ^ 



+ Hd w H(-i + cH)(H + H)(1 + ivH) 



(4.28) 



One may check that the Bianchi identities hold for either choice of sign for v. 

To compute g$ z , we first use the second equation in (3.3) to express g% z in terms of 
9iz, 92z and a and (3. The remaining a and (3 dependence can be expressed in terms of 
ip using (4.4). The result is 

3 1 1 

a 2 ± f3 2 /-[sinh(y9 — 0) + sinh(<p + 0) cosh(2(^) ± sinh 5 (2(p) sinh 5 (2ip)] 2 

o — = V2 3 1 (4.29) 

a P sinh*(2<p)sinh*(2<p) 
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Most of these quantities have already appeared in (4.24). The only new term is sinh((/?— <p) 
which is given by (4.6). Again one has to make choices for the signs of branch cuts. The 
intermediate steps are messy and we do not include them here. As before demanding 
the Bianchi identity to hold for g 3z restricts the possible branch choices. In this case 
however, one finds that the Bianchi identity for g 3z holds only for the branch with v — 1 
or equivalently Im(ff) < 0, except in the special case that c = 0. This is discussed further 
in section 5.2. The final result for g 3z is given by 



-1 r 



d w b 3 = -u 3 ih \ c\(—i + cH)(H + H)(l — HHf 



d w H(-i + cH)(l + H 2 ){H + H) 



+ d w H 8iH 2 H 3 + 3cHH 3 - iH 3 + 3cH 2 H 2 - \hiHH 2 - 2cH 2 + 2iH 2 H 



- cHH + 5iH + cH 2 - 3iH 



(4.30) 



4.3 Bianchi identities and equations of motion 

Several theorems exist guaranteeing that solutions of the BPS equations are automatically 
or almost automatically supergravity solutions, meaning that they are also solutions of 
the Bianchi identities and the equations of motion [3, 4]. In general, one needs to check 
the Bianchi identities and Maxwell equations, however the Einstein equations are then 
either automatic or automatic up to a single equation depending on whether the Killing 
spinors yield time-like or null Killing vectors. 

To make use of these theorems, we would first need to explicitly compute the Killing 
vectors and check whether they are null or time-like. Additionally, in deriving the explicit 
solutions for the metric factors and fluxes, we had to make sign choices for various branch 
cuts. These sign choices must be made consistently with the BPS equations, which 
requires checking if the BPS equations enforce additional constraints on the sign choices. 
The simpler alternative which we follow here is to check the equations of motion. We 
note that the equations of motion can still leave sign choices for the fluxes undetermined, 
with some choices breaking supersymmetry, and one must still check the BPS equations 
[11, 12]. 

In [1], the equations of motion and metric factors were checked for the special case 
c = and we use their methods here. We first make a change of notation and exchange 
the imaginary harmonic function h for a real harmonic function h with the identification 
h = —ih. Next, we note that the metric factors and field strength components define a 
solution only for G which satisfies its differential equation (3.24). In order to implement 
this constraint in general, we make a conformal transformation on £ to coordinates 
defined by h and its dual harmonic function, h. 



u = h + ih 
d u (h — ih) - 
2d u = d- h 



u = h — ih 
d a {h + ih) = 
2du = d h + id- h 



(4.31) 
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Equation (3.24) can be decomposed into its real and imaginary parts. 

G(h, h) =G R (h,~h)+iG I (h, h) 

d h G R + d h G r = ^ 
h 

d h G, - d~ h G R = (4.32) 

The second equation above demands that Gr and Gj are the derivatives of a single real 
function <f> and then the first equation is rewritten as a second order equation on <p. 

Gr = d h (<f>) Gi = d h (<f>) 

1(91 + 81-18^ = (4.33) 

Finally, we use this equation to eliminate all derivatives in h which are of order two or 
higher, which implements the differential constraint. To be specific, we eliminate <9|G, 
8fG and 8ld~ h G. 

h h n 

We begin with the Bianchi identity. 

dF = dvfiabi] - c.c. = i = {1,2, 3} (4.34) 

Using the expressions for d^bi given in (4.27), (4.28) and (4.30), one finds that the Binachi 
identities are automatically satisfied after first expressing G in terms of and using (4.33) 
to exclude all second and higher order h derivatives which appear. Next we check the 
equation of motion of the field strength. 

d* F + ^-F A F = (4.35) 

Using the anstaz given in section 2 and writing out in components we have 5 

= d^djdx + ^ I dwbidu, In I J -^- ) + c.c. ) f - I -f^— ) (d w b 2 du,b 3 - c.c.) 



2 



2 




3 



du,d w b 2 + ^ I du,b 2 d w In I J -^- ) + c.c. ) f - 1 -f^— ) (d w bidajb 3 - c.c.) 



3 



= du,d w b 3 + ^ ( du,b 3 d w In ( J -^- ) + c.c. ) - ^ ( -f^- ) (d^d^ - c.c.)(4.36) 

Again, one may check that these expressions are automatic after using (4.27), (4.28) and 
(4.30) for the fluxes and (4.19), (4.21) and (4.23) for the metric factors. We remind the 
reader that the Bianchi identity for g 3z selects the branch with i(H — H) > 0. We note 
that these equations are invariant under b± — > —bi along with a flip in sign of all of the 
metric factors so that /1/2/3 — > —71/2/3- Similar sign flips can be made for b 2 and 63, 
thus the Maxwell equations determine the fluxes up to an overall sign. We can determine 
the signs by considering the special cases c = {—3, 0, 3} where we take the geometry to 



5 We note that [1] contains a typo, where the factor of 1/4 in the third term of each equation should 
be a 1/2. 
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be either AdS 4 x S 7 (c = -3) or AdS 7 x S 4 (c = 0,3). This is discussed further in the 
next section. 

Finally, the Einstein equations can be checked. These are equations (9.19) and (9.20) 
of [1]. Some details of the derivation of these equations can be found there as well. 
The method for checking these equations is the same as with the Maxwell equations and 
Bianchi identities and requires replacing second and third order derivatives using 4.33. 



5 Summary of Solution and Remarks 

In this section we first give a summary of the general solution, gathering together all 
the formulas for the physical fields. We then give a discussion on the allowed range of 
the functions H and G. Finally, we show how AdS 4 x S 7 and AdS 7 x S* 4 fit into the 
framework of this paper and give examples of one parameter deformations as a function 
of c for each case. 



5.1 Summary of the general solution 

To help simplify the final expressions, we introduce the real harmonic function h related 
to the imaginary harmonic function h by h = ih. The solutions are then specified by the 
choice of a Reimann surface S, the real function h which is harmonic on S, the complex 
function G which satisfies the differential equation 

d w G= l -(G + G)d w \nh (5.1) 

and a choice for the value of c = (ci — c 2 )/c3, which can take any real value. The 
differential equation for G can be solved in terms of an integral after picking h as a local 
coordinate [1]. There is also a restriction on the range of allowed values for G. To express 
this constraint, it is convenient to use H — G/(l — icG), which, for any value of c, must 
satisfy the inequalities 

A\H\ A + (H - Hf < 

i(H-H)>0 (5.2) 

where the second inequality simply restricts the imaginary part of H to be negative. The 
range of allowed values for G is further discussed in section 5.2. The metric factors are 
given by 



p% _ \d w h\ ( 



c 4 h 4 \l + icH\ 



[1-\H\ 2 )(-4\H\ 4 -(H-H) 2 ^) (5.3) 



rr\ \ 3 



c 



3 h 2 (l-\H\ 2 ) 



(2\H\ 2 + i(H - H)) 
(4|#| 4 + (H - H) 2 y 



(5.4) 



.-7- \ \ 3 



^t2 



h 2 (l-\H\ 



(2\H\ 2 - i(H - H)) 



A\H\ A + (H 



Hf) 



(5.5) 
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The fluxes are given by 



c 3 



d w bi =2ui^zih 
c\ 



H(-i + cH)(H -3H + AHH 2 )d w H - HH(H - 3H + AHH 2 )d w H 
-ic({H - Hf - HH 3 + 3H 2 H 2 ^jd w H + H(-i + cH)(H + H){1 - iH)d w H 



(-i + cH)(H + H)[2HH — i(H — H 



(5.7) 



d w b 2 = - 2u 2 ^ih 
4 



H(-i + cH)(H -3H + AHH 2 )d w H + HH[H - 3H + AHH 2 )d w H 
ic({H - Hf - HH 3 + 3H 2 H 2 ^jd w H + H(-i + cH)(H + H){\ + iH)d w H 
-i ■ <U){U ■ ll)(2llll ■ nil I (5.8) 



4 



8iH 2 H 3 + 3cHH 3 - iH 3 + 3cH 2 H 2 - 15iHH 2 - 2cH 2 + 2iH 2 H - cHH 



+ 5iH + cH 2 - 3iH d w H + (-i + cH)(l + H 2 )(H + H)d w H + 



(—i + cH)(H + H)(l — HHf 



(5.9) 



where z/j = ±1. The equations of motion determine the fluxes up to an overall sign. In 
order to determine the signs, we must require the fluxes and metric factors to satisfy the 
BPS equations. This can be done for the special cases c = {—3, 0, 3} where we take the 
geometry to be either AdS 4 x S 7 (c = —3) or AdS 7 x S 4 (c = 0,3). Using the explicit 
formula for G and h given in section 5.3 to compute the fluxes and comparing to the ex- 
pressions given in section 3 of [1], we find v\ — 1, v 2 — — 1 and z/ 3 = —1. Additionally, as 
reviewed in section 2, the BPS equations admit a symmetry where one may flip the signs 
of any two fluxes. Conversely, when we flip the sign of a single flux the BPS equations 
do not admit any solutions. Thus we conclude that the configurations with {ui, u 2 , z/3} = 
{+, +, +}, {+, — , — }, { — , — , +}, { — , +, — } all preserve supersymmetry, while the config- 
urations with {ui, z/ 2 , z/3} = {+, +, — }, {+, — , +}, {— , +, +}, {— , — , — } all spontaneously 
break supersymmetry. This way of breaking supersymmetry is similar to the skew- whiffed 
solutions of [11, 12] 

5.2 Range of G and H 

It is interesting to determine the range of the functions G and H. Note that the range 
of H, depicted in Fig. 1, is fixed and independent of c (modulo a small technical point 
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-1.0 -0.5 0.0 0.5 1.0 

Figure 1: Range of H for any c. 



when c = 0, which will be explained shortly). As discussed above (4.14), the range of 
H is determined by the condition \H — H\ > 2\H\ 2 . The weaker condition \H\ < 1 is 
implied by this previous condition. In general, this allows for two branches of H, one with 
lm(H) > and one with lm(H) < 0. However, as discussed in section 4.2, the Bianchi 
identity for g 3z selects the branch with lm(H) < 0, except in the special case c = 0. 

The reason c = is special, is because the solutions corresponding to the mirror 
image range for which i(H — H) < are just an automorphism of the old range for which 
i(H — H) > 0. This can be seen as follows. An equivalent solution is one for which the 
solution G to its differential equation (3.24) is the same up to the symmetries discussed 
below (3.34) and the expressions for the physical fields are unchanged. Consider the 
mapping H — > —H which implies G — > —G and c — > — c. The change of sign in G 
is undetected by the differential equation which is also independent of c. However, for 
general c the expressions for the metric factors change in two ways: we must pick the 
other sign in front of i(H — H) and we must exchange C\ -B- C2 (where the latter is the 
only way to map c — > — c). For general c, this new combination differs from the original 
and is not a solution. In the expression for the physical fields, picking the other sign in 
front of i(H — H) undoes the negation of H but then it would be odd if switching c\ and 
C2 was a new solution in addition to the original one. However, c = means c\ = C2 and 
so the solution is simply mapped into itself. Thus the two branches are equivalent. 

We note that the restriction to v — 1 is consistent with the results of [1] . In the case of 
c = —3, corresponding to AdS^ x S 7 asymptotics, the range of G consisted of only a single 
branch. For the case of c = 0, corresponding to AdSj x S 4 , there were two branches for 
the range of G, one with Im(G) > and one with Im(G) < 0, which intersect at G — 0. 
One might wonder whether there are solutions where G interpolates between the two 
allowed regions, which intersect at G = 0. However, the explicit solutions constructed in 
[2] show that this is not possible, at least for the class of solutions considered there. 

The quantity H is convenient to work with precisely because its range is independent 
of c. However, the differential equation obeyed by H is not linear and depends on c. 
Since we will eventually have to impose boundary conditions for the differential equation 
of G, it is interesting to consider the range of G as well (which is most easily obtained 
from mapping the range of H). For c = 0, we have G = H and the range of G is of course 
simply the range of H given in Fig. 1. As c is taken to be more negative, the radius of 
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Figure 2: Range of G for c = —.25 (left), c = —.5 (middle) and c = — 1 (right) 
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Figure 3: Range of G for c = — 1.5 (left), c = — 3 (middle) and c = — 15 (right). 



the allowed range circle increases until it decompactifies at c = — 1 so that any G with 
a negative imaginary part is allowed. This sequence is illustrated in Fig. 2. Note that 
c = — 1 implies that c<i = which means that one of the spheres decompactifies. 

Continuing to decrease c from c = — 1, the range further increases such that all G are 
allowed except a disallowed circle whose radius decreases as c becomes more negative. 
This is shown in Fig. 3. Note that the middle figure in Fig. 3 is the range for c = —3 
which includes all solutions asymptotic to AdS^ x S 7 . Our range agrees with the one 
depicted in figure 2 of section 8 of [1] if we remember that our G is mapped into that of 
[1] by G c= __ 3 r]i = AG + i as discussed in (3.35). As c takes values in the range of {0, —3} 
the family of solutions interpolates, as a function of c between solutions asymptotic to 
AdSj x S 4 and those asymptotic to AdS^ x S 7 . As c — >■ —00 the range increases as seen in 
Fig. 3 (right) until it includes the entire plane. The case c = —00 corresponds to C3 = 
and is therefore the case of the other sphere decompactifying. 

We now consider what happens as c increases above 0. As c grows for c > 0, the 
range of G decreases as can be seen in Fig. 4. Note that the c = 3 case, show in Fig. 4 
(middle), contains solutions which are asymptotic to AdS? x S 4 but differs from the c = 
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Figure 4: Range of G for c = 1 (left), c = 3 (middle) and c = 15 (right). 



case in that the roles of the two 2-spheres are interchanged. Setting c = 1 implies C\ = 0, 
so this case corresponds to the decompactification of the AdS% region and is shown in 
Fig. 4 (left). 

In order to interpolate as a function of c between the two solution classes which 
asymptote to AdS 7 x S 4 corresponding to c = 3 or c = 0, one must pass through 
solutions with c = 1 for which the AdS^ part of the metric is decompactified. Similarly, 
to interpolate between solutions which are asymptotic to AdS^ x S 7 with c = —3 and 
solutions which are asymptotic to AdS-j x S 4 with c = 0, one must pass through the 
solutions with c = — 1 so that one of the spheres must decompactify. To interpolate 
between the AdS± x S 7 solutions with c = — 3 and the AdS 7 x S 4 solutions with c = 3, 
the other sphere must decompactify at c = oo. We note that this interpolation requires 
passing through the point c = oo and it is likely that one should study this region as a 
function of c -1 or some other well-defined parameter. There is no obvious relationship 
between the range of G and the full geometry of the solutions, but we expect that knowing 
this range will be necessary for finding the global solutions which we intend to work on 
in the future. 



5.3 AdS-j x S 4 , AdS± x ^S* 7 and a one parameter deformation 

Here we show how AdS^ x S 7 and AdS-j x S 4 fit into our general class of solutions and 
give one parameter deformations of both cases. We start with AdS 7 x S 4 , which is given 
by the following choice of parameters 

sinh(iL> — w) 
c = G = i- 



sinh(2w) 
Lc 2 

h = -i^cosh(2w) + c.c. H = G (5.10) 

where we have chosen the overall coefficient of h so that L is the the radius of the S* 4 . 
Writing w = x + iy, the domain of £ is given by < x < oo and < y < tt/2, with 
the boundary of AdS-j located at x = oo. One can check that H satisfies the constraint 
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(5.2). 6 

We can construct a one-parameter deformation of AdS-j x S* 4 as follows. Keeping h 
fixed, we first use the real scale and imaginary shift symmetry of the differential equation 
to write a more general G with two real parameters a and b. 

„ . sinhfw - w) 

G = ia , . _' + ib (5.11) 

smh(2w) 

Next we construct H for an arbitrary value of c using H = G/(l — icG) with the require- 
ment that H satisfies the same boundary conditions as the c = case so that H = when 
y = or y = n/2 and H = —i when x — 0. This leads to the restriction a = 1/(1 + c) and 
6 = with c left as a parameter. Note that for these choices, H satisfies the constraint 
(5.2). We leave a fuller discussion of boundary conditions to future work but below we see 
that the preceding simple choice of boundary conditions on H gives a regular geometry. 
The corresponding metric factors as a function of c are given by 

,6 _ r 2 1 6cosh 2 (a;)[l - c + (1 + c) cosh(2rr)][l + ccos(%) + (1 + c) cosh(2x)] 
Ji — L 



;i + c)(i- c 

_ inM - r - I I -- ,A r-ndif) <■)] 



2 16[1 -c + (1 + c)cosh(2x)][l + ccos(4i/) + (1 + c) cosh(2a;)] 



'1 + c) 4 cosh 4 ( 



X 



6 _ ^ 2 128cosh 2 (x)[l + ccos(%) + (1 + c) cosh(2a;)] sinh 6 (x) 

" ' (1 + c) 4 [l - c + (1 + c) cosh(2a;)] 2 

6 _ ^ 2 2cosh 2 (a;)[l -c + (1 + c) cosh(2a;)] sin 6 (2|/) 
^ ~ (1 + c)[l + ccos(%) + (1 + c) cosh(2x)] 2 



(5.12) 



One may readily check that the bulk geometry is regular when c is in the range — 1 < c < 
I. In particular, the S 3, s cap off smoothly at the boundary of S. When c is in the range 
— oo < c < —1, the geometry will contain a singularity, which can be seen by noting that 
fi will vanish at some value of x and y. The geometry is also smooth for c in the range 
1 < c < oo. The other AdS 7 x S* 4 solution with c = 3 is given by taking 

% sinh(u> — w) i 

c = 3 G = — — — — 

4 smh(2u?) 4 

t / d ~r 9x , \ t t sinh(u; — w) + sinh(2w) . , „. 

h = -i(2Lcl) cosh (2w) + c.c. H = -i — — ^— '- — ^ — 5.13 

v 3; v ; sinh(2w) -3sinh(w-w) y ' 

We now move onto the AdS± x S* 7 case, which is given by the following choice of 
parameters 

i cosh(w + w) i 

c = —3 G = — — — 1 — 

4 cosh(2w) 4 

i -Lc\ . , _ s TT . cosh(2w) — cosh(w + w) _ iX 

= ^sinh 2w + c.c. = % — - / / ^ 5.14 

4 v 7 cosh(2w;) + 3 cosh(w + w) 

where we have chosen the overall coefficient of h so that the radius of the S 7 is L. 
Writing w = x + iy, the domain of £ is given by — oo < x < oo and < y < ir/2, with 



6 We note that the G here is the negative of the G given in [1]. We remind the reader that in the 
special case c = 0, there are two equivalent branches of solutions, while only one branch extends to 
general values of c. 
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the boundary of AdS± located at x — ±00. In these coordinates the boundary of AdS± is 
split into two halves which are glued together along the AdS$ boundary. One can check 
that H satisfies the constraint (5.2). 

We can construct a one-parameter deformation of AdS& x S 7 using the same technique 
as before. Keeping h fixed, we first use the real scale and imaginary shift symmetry of 
the differential equation to write a more general G with two real parameters a and b. 

% cosh(u? + w) % . 1 . 

G = -a- y — — -)- — -)- %b (5.15) 

4 cosh(2w) 4 v ; 

Next we construct H for an arbitrary value of c using H = G/(l — icG) with the require- 
ment that H satisfies the same boundary conditions as the c = case so that H = 
when y — and H = —i when y = n/2. This leads to the restriction a = —2/(1 + c) and 
b = —(3 + c)/4(l + c) with c left as a parameter. Again, for this choice one can check 
that H satisfies the constraint (5.2). The corresponding metric factors as a function of c 
are given by 

ft = L 2 n 8C °^ X) , [ - 1 - c - (3 + c) cos(2 2/ ) + 2cosh(4:r) 
(1 -c) 6 (-l -c) L 

2 - (3 + c) cos(2y) - (1 + c) cosh(4a;) 

1 - c - (3 + c) cos(2y) + 2 cosh(4x)] [2 - (3 + c) cos(2y) - (1 + c) cosh(4x)] 



(2/?) L (l + c) 4 cosh 4 (2a;) 

6 _ 2 64 cosh 2 (2x) [2 - (3 + c) cos(2y) - (1 + c) cosh(4a;)] cos 6 (y) 
J2 — L 



[1 + c) 4 [l + c + (3 + c) cos(2y) - 2 cosh(4x)] 
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6 _ 2 8 cosh 2 (2x)[-l - c - (3 + c) cos(2y) + 2 cosh(4a:)] sin 6 (y) 
/3 (-l-c)[-2 + (3 + c)cos(2y) + (l + c)cosh(4x)] 2 1 ' ' 

One may readily check that the bulk geometry is regular when c is in the range — 00 < 
c < — 1. In particular, the S^'s cap off smoothly at the boundary of S. The geometry 
contains singularities when c is outside of this range, as can be seen by noting that f\ 
vanishes for some value of x and y. 

There are more solutions one can construct. For example using the same methods as 
above we can deform the solutions of [2] or the Janus solutions of [7]. It is also possible 
that new types of solutions exist, which carry both M2 and M5 brane charges when c is 
away from the special values c G {—3, 0, 3}. Another interesting possibility is to look for 
solutions which have AdS^ asymptotics. We leave the constructions of these solutions, as 
well as their implications for AdS/CFT to future work. 
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